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We study the magnetic properties of the 3d Hubbard model at half-filling in the TPSC formalism, 
previously developed for the 2d model. We focus on the Neel transition approached from the 
disordered side and on the paramagnetic phase. We find a very good quantitative agreement with 
Dynamical Mean-Field results for the isotropic 3d model. Calculations on finite size lattices also 
provide satisfactory comparisons with Monte Carlo results up to the intermediate coupling regime. 
We point out a qualitative difference between the isotropic 3d case, and the 2d or anisotropic 3d 
cases for the double occupation factor. Even for this local correlation function, 2d or anisotropic 3d 
cases are out of reach of DMF: this comes from the inability of DMF to account for antiferromagnetic 
fluctuations, which are crucial. 



I. INTRODUCTION 



The problem of strongly correlated fermions is a long-standing opened question. Recently progress in the resolution 
of the simplest model, namely the Hubbard model in rf = cx3 by the Dynamical Mean Field approach (DMF), 0] ||] 
marked a real breakthrough in the field. Nevertheless, the lack of spatial dependence and the difficulties in including 
the 1/d corrections, and the low dimensionality or the anisotropic character of the experimental systems, justify 
the ongoing research on more realistic dimensions d = 2 and 3. In low dimension, especially d = 2, one expects 
a strong k-dependence of the self-energy that is absent from d = oo calculations. In this paper we study the 3d 
isotropic Hubbard model, using a formalism previously developed in 2d, ||] [Q named TPSC for Two Particle Self 
Consistent. In 2d, the TPSC approach has been shown to give reliable results in comparisons with Monte Carlo 
calculations for one and two-particle properties, static and dynamical quantities, for a whole range of parameters: for 
a general band filling, on the whole temperature range accessible to Monte Carlo simulations (but probably not at 
very low temperature at half-filling in 2d), from low up to intermediate value of U/W, for the nearest neighbor model 
as well as for the Utt' case. Q The strengths of this approach rely on the respect of several benchmark properties: 
it fulfills the Pauli principle and is conserving in the Baym and Kadanoff sense. |^ It was shown in 2d to satisfy the 
Mermin- Wagner theorem. Furthermore the TPSC approach is appealing due to its simplicity: the vertices retained 
in the particle-hole channels are static and local. 

Apart from the 2d case, the highly anisotropic 3d case (inter-plane hopping much less than in-plane one) has been 
studied previously. Q It was found that the antiferromagnetic order can be established in the regime of incoherent 
3d one-particle motion and the crossover between 2 and 3d was studied in details. In this paper we want to explore 
the magnetic properties of the isotropic 3d model, focusing on the paramagnetic phase and on the Neel transition 
approached from the disordered side. We make comparisons with other approaches, including Monte Carlo results, and 
DMF approach for the hypercubic lattice. We also stress the qualitative difference between isotropic and anisotropic 
3d cases, even for a local function as the double occupation factor. 

The paper is organized as follows: In Sec. II we review the TPSC formalism for the particle-hole properties, namely 
spin and charge susceptibilities. Sec. Ill concerns the Neel transition, and Sec. IV is devoted to the temperature 
dependence of the double occupation factor from which all the magnetic properties follow. Finally we summarize our 
conclusions in Sec. V. An appendix details the calculation of the double occupation factor for a simple two-site model 
that enlightens the high-temperature lattice results. 

II. TWO-PARTICLE SELF-CONSISTENT FORMALISM 

We start from the Hubbard Hamiltonian, 

H = - ^3 {4<TCja +c]^Ci^^ +Uj2^i'i^il ' l^J2 ^^'^ ' 

<ij>(T i i<T 
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where the operator Ci^ (cj^.) destroys (creates) an electron of spin a at site i, rii^ is the density operator, and tij is the 
symmetric hopping matrix. The on-site repulsion U is the screened Coulomb interaction, fi the chemical potential. 
For a while hopping parameters will be kept general. Indeed the Two-Particle Self-Consistent formalism (TPSC) 
approach only requires a local interaction but can be applied to a general band structure. 

Let us sketch the TPSC approach, fully developed in Ref. ||l and |]. The TPSC can be formulated in a closed 
form in the functional integral formalism. Here we shall concentrate on the two-particle properties, and obtain the 
useful relations in the equation of motion approach. In the presence of a time, space and spin varying external field 
(j), a generalized Green function can be defined: B 



Tr 



e~f^"TrSiP)cUn)c,A^ 



G{n,r,,n;q^)- —^-^^^—^^ , (2) 

where means imaginary time ordering and 

ln^(/3) = -/VVct^(r)c,,(r)0,(/,r) . (3) 

All the field dependence is in S{P) and 

cl{r)^e^"cle-^'' ■ (4) 

To obtain the spin and charge susceptibilities, the equation of motion for G has to be evaluated. Then in the linear 
response theory, the functional derivative with respect to <p is taken at zero field leading to: 



with the kinetic part 



kq 

All the operators in Eqs. (^) and are at the same time ti and (...) is a short-hand notation for averaging with 
S{I3) as in Eq. (H). ct is opposite to ct, k is a wave vector of an iV-site lattice, and is the Fourier transform of t^j-. 

The ansatz of the TPSC approach is to write, for the term proportional to U in Eq.(||), that the following approx- 
imation retains the main correlation effects: 

S (^^a^jcr'^ls') — ^ (^9tTail,Ti) l^cl^Cja'^ (nia)^ for 1^1 Or l^j , (7) 

with the following definition: 

gaaihn)^- — — — - (8) 

{nia(Ti)) {nwiTi}) 

Eq. ^ is exact when i — j- The expansion of the functional derivative with respect to cj) leads to 



since at zero field in the paramagnetic phase, the coefficient of the term proportional to 5 ycJ^Cjo-y cancels. 

To get the spin susceptibility it remains to take (paii, Ti) — <J(i>{i, t;), then to time- and site-Fourier transform Eq.(^, 
using Eqs.(0) to (0). From the definition 
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the preceding equations lead to an effective RPA expression for the spin susceptibihty 



with 

Vsp^gnV ■ (12) 

gil is the time- and site-independent version of Eq.(^). This q and iujq independence foUows from the cancellation 
of the term proportional to Sgcrg which does not depend on the spin index, and disappears by summation on a: 
X^cr ^^9o-a = 0. As will be recalled later, g^i can be evaluated in a self-consistent equation dictated by sum rules. 
For the charge susceptibilities the term proportional to Sg^a doesn't cancel, and one obtains a q and iwg-dependent 
renormalized interaction. However further simplification consisting in taking a q and iwg-constant renormalized vertex 
for charge Uch = {9u + ^g)U , where 5g can be self-consistently evaluated using exact sum rules, has been shown to 
give reasonable results in comparison to Monte Carlo simulations. The sum rules to determine self-consistently Ugp 
and Uch follow from the the Pauli principle and the fluctuation-dissipation theorem, here expressed in Matsubara 
formalism 

{St St) = {[n, - n,Y) =n-2 {n,n,) = E . i^;/^'7\ . , (13) 



and 

I3N 1 + 5C/chXo(q, iUJq) 

where (3 = 1/T, n = (n|) -I- (n^), and Srii = rii^ — n. The necessity to have two different irreducible vertices for 
spin and charge is imposed by the Pauli principle. Simple RPA for which the spin and charge susceptibilities differ 
only by a sign in the denominator cannot fulfill these previous rules. 

Although the susceptibilities have an RPA functional form, the physical properties of the theory are very different 
from RPA because of the self-consistency conditions for determining Ugp and Uch . This approach has been shown j|] 
to take into account both the local quantum rcnormalization effects (Kanamori-Brueckner screening) which reduce 
the Coulomb repulsion, and long-wavelength thermal fluctuations, giving a phase transition only at zero-temperature 
in two dimensions. This theory is self-consistent at the two-particle level (hence its name), in contrast with FLEX |^ 
for example, which is self-consistent at the one-particle level. 

In this paper we shall concentrate on magnetic properties for the nearest-neighbor model on a 3d-simple cubic 
lattice with a spacing a — I, for which the dispersion relation reads 

3 

ek = — 2t E cosfci, (15) 

1=1 

with a bandwidth W = 12t. The numerics only concern the half- filled band case, hence Xo(cIjO) is maximum at 
Q = (tt, tt, tt) due to the perfect nesting. The calculations are done on a variable mesh for the Brillouin zone: 16^ on 
the coarse one, 16'^ more points near Q, but as will be detailed below, some integrals are analytically evaluated, and 
the non-interacting susceptibility is finely evaluated at low temperature. We also present calculations on finite-size 
lattice for comparisons with Monte Carlo results. Hereafter t will be the energy unit, and we set fcs = 1. 



III. THE ANTIFERROMAGNETIC TRANSITION 



The TPSC formalism has been shown to give reliable results far from the transition and also in the part of the critical 
regime that is accessible to Monte Carlo calculations. Important constraints such as the Mermin- Wagner theorem 
arc fulfilled by the TPSC approach. Indeed in 2d as the temperature goes down, Ugp decreases, thus preventing a 
finite-temperature magnetic transition. However deep in the low temperature regime at half-filling in 2d, the relation 
Usp = U g^i breaks down: there is no reason to believe that there is no double-occupied site in the ground state. 0| In 
the highly anisotropic three-dimensional case with a small inter-plane hopping tz <ti t, a finite-temperature transition 
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can be stabilized. Here we shall concentrate on the isotropic case and make comparisons with other approaches 
developed for the 3d or for the infinite-dimensional Hubbard model. 

In 3d, at half filling, the self-consistent Eq.(|l3[) from which the magnetic properties will be evaluated, can be 
rewritten 



(16) 



where Xspi'ij 0) is the asymptotic form of Xsp(q, 0) obtained from a Taylor expansion around Q of the denominator 
and defined on a domain of radius A around Q. This enables an analytic evaluation of the diverging part of the 
susceptibility. Defining Umfc = 2/xo(Qj 0) the mean-field critical value of U at which the transition occurs in a RPA 
treatment, and 6U = Umfc — Usp the small energy scale that measures the proximity to the phase transition in the 
TPSC approach, the asymptotic form of the spin susceptibility simply reads 



1 



5U 



i+e 



where ^ is the magnetic correlation length defined in terms of a microscopic length S.o' 

f2 _ c2 ^sp 



(17) 



(18) 



d'xo 



dql 



(19) 

q=Q 



The Neel transition is reached from above when Usp — > f^m/c- ^^'^ TPSC approach validity is restricted from low up 
to intermediate value of the ratio U /W . Indeed it has been shown, and will be confirmed below, that Usp saturates 
at strong coupling while one would expect Usp to have a maximum as a function of [/, as suggested by RPA at low 
repulsion and by the Heisenberg model limit at strong repulsion. Nevertheless the reliability for U ^ W va 2A has 
been established by comparisons with Monte Carlo calculations, and gives some confidence in 3d for the intermediate 
coupling regime. 

Fig.l presents Tjv as a function of U in units of the hopping parameter t = \. TPSC results are presented along with 
Monte Carlo, RPA, dynamical mean field results on an hypercubic lattice, |11| ||l2| spin-fluctuation formalism. 



I3| and spin Hamiltonian approach. ||1J] ||15[] The last one is a mapping of the Hubbard model onto an Heisenberg 
model with arbitrary-range, U-dependent exchange integrals Jij{U), treated in an Onsager reaction field method. The 
DMF results are scaled for comparisons with id by ^< e^, > (details will be given in next section). 
Usual RPA predicts an exponential behavior: 



where P = 3.8 ± 0.1, and iV(0) ~ 0.14 is the 3d DOS at the Fermi level. Due to the same functional form, the TPSC 
formalism also predicts an exponential form for Tm with U replaced by Usp- The saturation for [/ > of Usp to the 
value ~ 3.5 leads to the saturation of T/v as can be seen in Fig.l. Usp = Ug^i is not the correct ansatz at large U. 

Recently there have been a few controversies concerning the Neel temperature. One of these concerns the relevance of 
the RPA at weak coupling. |16| ||l^ A discrepancy was found in an infinite-dimensional approach with 1/d corrections 



in Ref. [16|, and in pure 3d formalism in Rcf. [|7|. From Fig.l it appears that the TPSC curve merges with RPA 
results, but only at very small coupling. The modified Stoner criterion used to obtain the Neel temperature in Rcf. 
Ip^l uses a maximally crossed diagram series which was shown to be not reliable close to half-filling in 2d. jl^ Up to 
U ~ W, TPSC and DMF approaches are very close. 

The second controversy concerns the interpretation of the Monte Carlo results: Hasegawa [|l^ has argued that 
the Neel temperature extracted from the Monte Carlo simulations was overestimated, in fact it would be rather a 
mean- field temperature than a true T/v. A indication of this can be seen at large U, where the Monte Carlo is closer 
to the mean-field result T¥^ ~ /U , than to the high temperature series result T^-^ ~ 3.83i^/C/. However, as 
pointed out by Hirsch pC| this could be also an indication that the itinerant nature of the Hubbard model persists 
up to large coupling, enhancing the transition temperature. We cannot settle this debate at large U , the TPSC being 
inappropriate in this regime, but even at low interaction, our results are not very close to the Monte Carlo points. 
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The U dependence of T/v can be traced back using the self-consistency equation (^^. Close to the phase transition 
the most important contribution is taken into account by the asymptotic form of the susceptibility. At T/y 5Ti ~ 
Umfc/U, and xip simply reads 

X%{l.^)-T^^ , (21) 



for q measured from Q. It then follows that the and U dependence of Eq.(|l6|) can be made explicit 

l^-R^T^j^jj^A , (22) 

where A is the radius of the sphere centered on Q where Xsp is defined, and i? is a short-hand notation for all terms 
which weakly depend on T/y and U. It has been shown that at low temperature °c scales as T^^, in the 
two-dimensional case. This also holds for the 3d case where 



d^Xo 



c 

Q ^ 

with c ~ 0.0705. This behavior extends over a wide temperature domain including T/y. It then follows that 



(23) 



. (24) 

On the whole range of T/v, the scaling behavior of U,nfc — 2/xo(QiO) cannot be neglected and is set, for T < t, by 
Xo(Qj 0) — 2A^(0) \n{t/T) + (3 where [3 = 0.38. Equation ( p4|) establishes the Neel temperature saturation. 

The TPSC formalism has been shown in 2 < d < 4 dimensions to be in the same universality class as the Berlin-Kac 
spherical model, i.e. in the same class as the Heisenberg N-component model with N oo. The corresponding 
critical exponents in o? = 3 are: v = 1, j = 2, rj = 0, a = —1, P = 1/2, 6 = 5. Furthermore the dynamical exponent 
has been shown to be z = 2. These exponents were also found in the Spin-Hamiltonian approach treated in an 
Onsager reaction field theory. j2^ The size of the critical region, i.e. for which ^ ^ 1, is quite narrow. Well above 
the critical regime, the behavior of with T remains roughly linear on the whole temperature range studied up to 
T ^ 2t, even for small U, and does not show any crossover to the RPA regime characterized hy i/ = 1/2. |2^ The 



linear behavior of ^ ^ at high temperature can be ascribed to (Eqs.(|lq), ([191 ) and (23)). However for high T, the 
temperature behavior of xIpiQ' 0) is the mean-field one. Indeed, from Eq. (|ll[), we have 

x;p{Q,o)^xo\Q,o)-^Ugn ■ (25) 

At high temperature g-^i slowly increases with T (see next section), and x7p essentially behaves as Xo^- 

To conclude this section we show in Fig. 2 the value of the local moment /i^ (5,f 5f) at Tjv as a function of U. 
TPSC results are compared to DMF results. [|ll[ ||l^ Due to the saturation of Ugp, the moment does not fully develop 
at large U in contrast with the DMF results and with what is expected from the Heisenberg picture. Nevertheless, 
up to C/ ~ 8, TPSC gives a very good agreement compared with d = oo results. 



IV. PARAMAGNETIC PHASE PROPERTIES 



The central quantity that determines the magnetic properties of the TPSC approach is the local correlation factor: 
(7ll = (n^ni) I (ri|) (rij), which is simply related to the fraction of doubly occupied sites and to the local magnetic 
moment {S'f) = 3(n — ^n^g-^i). Furthermore, g^i is closely related to the self energy through the following identity 

i 

^Trl:G=U{n^n^) , (26) 

where G and S are respectively the one-particle Green's function and self-energy. From this exact relation, and in 
spite of the local nature of gi^i, one expects a strong dimensional effect due to the influence of the magnetic fluctuations 
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on the self-energy S in 2d. Q] We are primarily interested in the isotropic 3d case, but we will also comment upon 
the 2 and anisotropic 3d cases. 

In Fig. 3 the temperature dependence of g-^^ predicted by TPSC for 3d, is displayed for various values of U. As 
expected, increasing U is very effective to reduce the fraction of doubly occupied sites. For U W, Usp saturates 
and thus g-fj^ ^ \/U while it should decrease more rapidly with U. In this regime the TPSC approach breaks down. 
At high temperature g-^<^ slowly increases towards 1, whereas on the low temperature side, we are limited to T > Tjv. 
The arrows locate Tmin, the minimum of the curves. 

g^l has been evaluated in the DMF approach. [Q, |^ To make quantitative comparisons between TPSC and DMF 
results, one needs to use the kinetic energy scale y/< e''j[ >. Indeed, within the infinite-dimensional model on the 
hypercubic lattice, the hopping amplitude is scaled to keep U j ^ < ef, > finite, whereas the bandwidth is infinite. 

For a finite dimension d, averaging over the Brillouin zone leads to \/ < > = \/2d- In the comparisons of the 
coupling dependence of g^ i at two temperatures, a quantitative agreement is found between the TPSC approach and 
the DMF results as reported in Fig. 4. The TPSC curves correspond to T = 0.6 and T = 1.7, i.e. they surround T„ 



whereas the infinite-dimensional results extracted from Ref. |2^, where < e| >= 1/2, are for T = 0.175-\/l2 ~ 0.6, 



^ ^ - I/O „„„ Ti _ n 1 -TK. /To ^ 

and T = 0.5\/T2 ~ 1.73 respectively. Up to half the bandwidth {Ud=oo ~ 2 in d = cx) units), the numerical values are 
quite close. 

In both DMF and 3d-TPSC approaches, a non-monotonous temperature behavior is observed: for U not too strong, 
gijl first begins to decrease before increasing at higher T. The value of the temperature corresponding to this minimum 
Tmin marked by an arrow in Fig. 3, is plotted in Fig. 5 and compared to the infinite-dimensional results. ]2^ ], |^ They 
are quite compatible, since the error bars in Fig. 5 for the TPSC results come from the flatness of the curves. For 
DMF, Tmin ^ for J7 ~ 4-\/T2 ~ 14, whereas in 3d for U ^ 10, gn increases monotonously for T > T/v, and cannot 
be calculated within our approach for T < T^. As discussed by Georges et Krauth, |Q the decrease of {n-^n,]) reflects 
an incipient localization tendency when the temperature increases. This type of behavior is made more acute in the 
Mott transition, observed for example for the V2O3 compound doped with Cr: heating can drive the transition from 
metal to insulator ||2^ (but the hypothesis of a structural effect was raised 26 1). The same idea is at work in the 
liquid-solid transition in ^He, in the so-called Pomeranchuck effect: the spin entropy contribution is greater in the 
solid phase, and increasing the temperature can drive a transition to the solid state. 

From Eq.(p6|) a strong dimensional effect is expected for g^j^, and it becomes interesting to compare the TPSC 
predictions for 2, 3 and anisotropic 3d cases to the D MF results. To make comparisons between different flnite 
dimensions d, we used again the energy scale \J < ej.> rather than naively W ^ (we flrst used W and obtained a poor 
agreement). yj< e\> is a better measure of the kinetic energy: it takes into account not only the bandwidth, but 
also the particular shape of the density of states. Furthermore, this is the same scaling as the one previously used 
for comparisons between d = 3 and d = 00. In the anisotropic 3d case, the small inter- plane hopping <^ t does not 
strongly affect the kinetic energy, and as for the 2d case one has •\/< e^. > = 2. 

In Fig. 6, we show the temperature dependence of gi^i for the two, isotropic and anisotropic three-dimensional cases 
in the TPSC approach. The units on the figure are the 3d ones. We flrst focus on the 2 and isotropic 3d cases. 
The values of gi^i cannot be distinguished for 2 and 3d at high temperature. At a temperature close to the arrows 
locating the minimum of the isotropic 3d case, the curves separate. At low temperature in 2d, Ugp goes down when T 
decreases. This behavior is conflrmed by Monte Carlo calculations for the double occupation factor, |2^ and is a 
manifestation of the incipient antiferromagnetic fluctuations, as can be seen by inspection of Eq.(13): the increase of 
the spin susceptibility integral, leads to a decrease of the double occupation factor. This behavior is not clearly seen 
in the 3d case due to the very narrow range of the critical regime, and to the smaller effect of fluctuations in high 
dimensions. The temperature at which the 2 and 3d behaviors separate, close to Tmin, is higher than the crossover 
temperature which marks the entering of the critical regime characterized by a large correlation length, p9| i.e. ^ does 
not need to be large to see a strong dimensional effect. The non-monotonous temperature behavior of gi^i observed 
in 3d and d = 00 cases, is also present but much less important in the 2d case for small interaction values, as revealed 
by a careful inspection of Fig. 6 for U = 2. In 2d the strong antiferromagnetic fluctuations prevent a sensible increase 
of g-\i when the temperature decreases. 

In the anisotropic 3d case, a small inter-plane hopping can stabilize the ordered phase, but has no effect on the 
value of gii, as shown in Fig. 6 for ?7 = 4, and tz = 10~^ (3d units): for T > T^itz) ~ 0.18 the 2d and anisotropic 3d 
curves cannot be distinguished. This behavior is a manifestation of the fact that for these U and tz values, the Neel 
transition takes place in the incoherent 3d one-particle motion regime: the first fermionic Matsubara frequency ttT/v is 
much higher than tz- The inter-plane hopping having no effect on the one-particle Green's function and self-energy, it 
follows from Eq. (p6|), that the behavior of is dominated by two-dimensional effects, even if the critical exponents 
are the same for the isotropic and anisotropic 3d models. 

In the high temperature regime (typically for T > Tmin), when the dimension does not have any influence, one 
can expect that a simple two-site problem may be instructive. The calculation of gn in the two-site case is detailed 
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in the appendix, we list here the main results. The two-site problem solved in the Grand-canonical ensemble, with 
a chemical potential fixed to the half-filling value, gives a good quantitative agreement for the double occupation 
factor in comparisons with the TPSC, at high temperature, for low and moderate interaction values. The two site 
results also display a non monotonous temperature behavior for g^i ^ with a {/-dependent minimum. This minimum is 
pushed to high temperature when U increases, in contrast with what was observed in the TPSC and DMF approaches 
(Fig. 5). However if we freeze the charge fluctuations which arc important in a 2 site problem, working in the canonical 
ensemble, the decrease of Tmm when U increases is recovered. 

From the comparisons between 2d, 3d, d = oo, and the two site case, it appears that the non- monotonous temper- 
ature behavior of (n|n|) is not a purely dimensional effect, nor a strong U effect, and that it is overcome by strong 
antiferromagnetic fluctuations in d = 2. 

Finally, we compare our results on finite-size lattices with Monte Carlo calculations. In Fig. 7, we compare our 
results for fi"^ = (SfSf) to the Monte Carlo results of Hirsch [|o) on a 4^ lattice at T = 0.5. The TPSC calculations 
are also done on a 4^ lattice. Clearly up to half the bandwidth, the results are quite compatible. In Fig. 8 we have 
plotted the staggered susceptibility as a function of the site number (4,6,8 and 10^), for different temperatures and 
for [/ = 6. These are compared with the Monte Carlo results of Scalettar et ai. An overall agreement can be 
observed except at T = 0.4 for the intermediate sizes. In the TPSC and in Monte Carlo calculations, a size effect is 
clearly seen at this temperature, however for Monte Carlo, the linear dependence of Xsp(QjO) is the signature of an 
"ordered" phase, whereas the transition is not reached yet in TPSC. 

V. CONCLUSION 

The TPSC approach has been shown in the isotropic 3d case, to give compatible results in comparison with d — oo 
calculations, for the Neel temperature as well as for the local correlation function (n^ni). The comparisons with 
Monte Carlo calculations from low up to intermediate coupling regime, except for the Neel temperature which is 
systematically larger in Monte Carlo, are remarkable. 

The various results for obtained primarily for 3d, but also compared to 2 and anisotropic 3d Hubbard models 

shed light on the applicability of the DMF approach. The latter compares very well with the isotropic 3d case for 
(n|n|), which is a local quantity that depends on the dimensionality. That may be seen through Eq.(p6|). However, 
due to the strong k-dependence of S, the 2d case, or the anisotropic 3d case for which the Neel transition can be 
reached with a purely two-dimensional self-energy, are out of reach of the d — oo model, as long as the l/d corrections 
are not taken into account in a satisfactory manner. Progress in that direction has been achieved recently. ]30|] 

Finally, let us stress the fact that up to intermediate coupling, the TPSC gives the same agreement in 2 and 3d, 
and thus is able to capture the essence of two very different physical systems. This method is not restricted to 
half-filling, nor to the hypercubic case, in contrast with many other approaches. It then follows that the frustrated 
antifcrromagnetism as well as the metallic fcrromagnetism can be considered. Furthermore a generalization of this 
approach can be elaborated in an ordered phase. Work in that direction is in progress. 

It's a pleasure to thank A.-M. S. Tremblay, J.-L. Richard and S. Pairault for enlightening discussions and key ideas. 
We are indebted to Yuri Vilk for providing us with his code enabling the calculations in 2d, and to Steve Allen for 
sharing his Monte Carlo calculation results. 

APPENDIX A: TWO-SITE RESULTS 

In this appendix we detail the calculations of the double occupation factor for the 2 site problem. From the compar- 
isons between the 2d, isotropic and anisotropic 3d, and d = oo results, presented in the main text, and the observation 
that, when properly scaled, (n-^ni) does not depend on the dimension at high T, we expect a good quantitative agree- 
ment between the 2 site problem and the infinite lattice results for any dimensions at high temperature. Apart from 
a quantitative agreement, it will be shown that this simple problem roughly explains the non-monotonous behavior 
of gii as a function of T, and the [/-dependence of Tmm observed in the 3d and d = oo lattice case. 

Considering a two-site problem, with a chemical potential ^ = U2 sites I fixed to have a mean occupation number 
per site n = 1, it is a straightforward exercise to diagonalize the Hamiltonian, and to evaluate the double occupation 
factor. The results for g-\^i for this problem are presented in Fig. 9 as a function of temperature and various values of 
the interaction, along with the 3d TPSC results. As previously done for comparisons between different dimensions, 
interaction and temperature can be rescaled for quantitative comparisons. With the periodic boundary conditions 
for the two sites problem, the kinetic measure is \/ (ef.) — 2, as in 2d. For low and moderate repulsion, and for a 
temperature typically higher than the 3d Tmim both results are close. The 2 site case not only provides a rough 
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quantitative agreement, but also displays a non-monotonous behavior as a function of temperature. This one can be 
easily understood in terms of eigenvalues e and eigenstates. In Fig. 10 we have plotted the different values of e — ^.N 
as a function of C/2 sites- Table 1 gives a schematic description of the corresponding states and their contribution to 
("T"i)- 



Index 


degeneracy 




{(t)\n^ni\(j>) 


a 


1 


^ - +64i2 - 2^ 


1 1 u 


4 4 ^U'^+eit^ 


b 


3 


-2/x 





c 


2 
2 


-2t + U - 2,fi 

-2t- fi 


1/2 



d 


2 
2 


2t + U -3^J. 
2t- 1.1 


1/2 



e 


1 
1 
1 


U-2fi 

2U~4:fl 




1/2 
1 



f 


1 


^ + iVC/'+64i2-2Ai 


1 ; 1 U 


4^4 ^U'^+eit'-' 



TABLE 1. Details of Fig. 10. Index refers to the different letters indexing the curves in Fig. 10. For each curve we 
specify the degeneracy order, and the value of e — /iiV, the number of electrons being the coefficient of fi. The last 
column gives the value of the double occupation operator in each state. For example the curve c corresponds to 2 
three-electron states, with ((/)|n|nJ0) = 1/2 for each, and to 2 one-electron states with no double occupation. 

For U2 sites < 4, the first excited state is four-fold degenerate and corresponds to charge fluctuations (the particle 
number is not the same as in the ground state). Two of these four states correspond to a non-zero double occupation. 
For higher interaction the first excited state is three-fold degenerate, corresponding to spin fluctuations, and has 
no double occupation. Thus for U2 sites > 4 the simple 2 site problem predicts a non-monotonous temperature 
dependence of (n-^ni). In that case the distance between the first and second excited sates raises with interaction as 
can be seen in Fig. 10. It follows that the temperature at which the double occupation factor begins to increase, grows 
with the repulsion. This contrasts with the behavior observed in d = 3 or d = cx). Nevertheless, the simple 2 sites 
problem favors charge fluctuations compared with the thermodynamic limit, and interestingly the decrease of Tmj„ 
with an increase of U, is recovered if we freeze charge fluctuations taking 2 sites with precisely 2 electrons. If freezing 
charge fluctuations gives a qualitative agreement for T„i„ as a function of U, the quantitative agreement for g-^^ in 
comparisons with the infinite lattice case is not as good as it was with charge fluctuations allowed. 



^ For correspondence, email: dare@lrc.univ-mrs.fr 
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FIG. 1. Neel temperature as a function of U , for the TPSC formalism compared with other approaches, see text. 
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FIG. 2. = (S!S!) at Tn as a function of U. Black circles correspond to the d = 3 case, while open symbols are DMF 
results, extracted from Ref. [11]. The energy unit is the 3d one. 
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FIG. 3. Double occupation factor g-^ = (n|n|) / (n|) (n^) as a function of temperature for different values of the on-site 
repulsion. The arrows locate the minimum of the curves. 
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FIG. 4. Comparisons between 3d-TPSC and DMF for as a function of U for two different temperatures. Black circles 
correspond to the d = 3 case, while open symbols are DMF results, extracted from Ref.[23]. The energy unit is the 3d one. 
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FIG. 5. Temperature corresponding to the minimum of flti(^) the 3d Hubbard model compared to the DMF results 
extracted from Ref.[23]. Black circles are for 3d-TPSC, open symbols for DMF results. The energy unit is the 3d one. 
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FIG. 6. Comparisons between d = 3, d = 2, and the anisotropic 3d cases for (Jn as a function of T, and different U values. 
The energy unit is the 3d one (see text). Solid line corresponds to 3d, the dashed one to 2d. Square symbols are for the 



anisotropic 3d case, with tz 
curves. 



10 and Z7 = 4, this curve ends up at TN{tz) ~ 0.18. The arrows locate the minimum of the 3d 
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FIG. 7. fj? = (SfSi) at T = 0.5 as a function of U. Black circles correspond to the TPSC results on a 4^ lattice, open 
symbols are the Monte Carlo results on a 4^ lattice, extracted from Ref . [20] 
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FIG. 8. Xsp (Q:0) at = 6 and different T values, as a function of the number of lattice sites. Symbols with error bars, 
linked by dashed lines are Monte Carlo results from Ref.[10], while symbols linked by solid lines are the TPSC results. 
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FIG. 9. g-ji as a function of temperature for different values of the on-site repulsion (solid line) compared to the two-site 
result (dashed line). 
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FIG. 10. Eigenvalues of the two-site problem e 
be found in Table 1. 



nN as a function of the on-site repulsion. Details about these states can 
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